Striped States in Quantum Hall Effect: 
Deriving a Low Energy Theory from Hartree-Fock 



Anna Lopatnikova' 1,2 ', Steven H. Simon' 1 ', Bertrand I. Halperin' 2 ', and Xiao-Gang Wen' 3 ' 
™ Lucent Technologies Bell Labs, Murray Hill, New Jersey 07974 
' 2 ' Department of Physics, Harvard University, Cambridge, Massachusetts 02138 
' 3 ' Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 

There is growing experimental and theoretical evidence that very clean two dimensional electron 
systems form unidirectional charge density waves (UCDW) or "striped" states at low temperatures 
and at Landau level filling fractions of the form v = M + x with 4 < M < 10 an integer and 0.4 < 
x < 0.6. Following previous work, we model the striped state using a Hartree Fock approach. We 
construct the low energy excitations of the system by making smooth deformations of the stripe edges 
analogous to the construction of edge state excitations of quantum Hall droplets. These low energy 
excitations are described as a coupled Luttinger liquid theory, as discussed previously by MacDonald 
and Fisher (Phys. Rev. B 61, 5724 (2000)). Here, we extend that work and explicitly derive all of 
the parameters of this low energy theory using a Hartree Fock approach. We also make contact with 
the equivalent low energy hydrodynamic approach of Fogler and Vinokur (Phys. Rev. Lett. 84, 5828 
(2000)) and similarly derive the parameters of this theory. As examples of the use of these results, we 
explicitly calculate the low-energy excitation spectrum and study tunneling into the striped state. 



I. INTRODUCTION 

Even before the discovery of the quantum Hall effect, 
it was suggested that the ground state of a partially filled 
Landau level of a two dimensional electron gas (2DEGi) 
in a high magnetic field should be a charge density waveEl. 
It was clear that (neglecting disorder) due to the degener- 
acy of the noninteracting electron spectrum, the ground 
state of such a partially filled Landau level should be de- 
termined entirely by electron-electron interactions, and 
charge density wave states were obvious candidates. In 
fact, when the fractional quantum Hall effect was discov- 
ered for the one third filled Landau level, it was initially 
identified as a charge-density wave phenomenon. How- 
ever, it was quickly determined that such observed frac- 
tional quantum Hall states in the lowest Landau level are 
uniform density fluids. (States that appear to be Wigner 
crystals were later observed in the lowest Landau level at 
still higher fieldsB). 

In partially filled higher Landau levels, however—there 
is now growing experimentalota and theoreticalETE£l evi- 
dence that charge density wave states do occur at low 
temperatures in ultra-clean samples. The experimental 
phenomenology of these high Landau level states, which 
has been, elucidated by a series of publications by two 
groupsou, is quite rich and complex. For filling fractions 
v = nifio/B (with n the electron density, B the mag- 
netic field, and 0o = hc/e the flux quantum) between 
roughly 9/2 and 21/2, certain trends have been observed 
that may indeed be associated with charge density wave 
behavior. In particular, for filling fractions of the form 
v = M + x with M an integer between roughly 4 and 10 
and 0.4 < x < 0.6 an anisotropic state is observed that 
has been tentatively identified as a unidirectional charge 
density wave (UDCW) or "striped" state. 

Although the discovery of this state attracted quite a 



bit of attention in the community, perhaps it should not 
strike us as surprising. Indeed, such a state had been pra* 
dieted for the higher Landau levels seyeral years earliem 
by Koulakov_ Fogler, and ShklovskiiQ and by Moessner 
and ChalkerQ. Koulakov et al. found that, within the 
Hartree-Fock approximation, when the partially filled 
Landau level is nearly half-filled, the ground state of 
a 2DEG is a UCDW that looks like a set of stripes in 
momentum space. The filling factor in the stripes alter- 
nates between [v] and \v] + 1 ([X] is the integer part of 
X) with a period, a, about 2.7 times the cyclotron ra- 
dius R c (R c = \p2Jj + 11, where L is the partially filled 
Landau level, L = [M/2], and I is the magnetic length, 
I = (Hc/eB) 1 / 2 ). Moessner and Chalker confirmed the 
result for finite temperatures and justified the use of the 
Hartree-Fock approximation by proving it to be exact in 
the limit of infinite Landau level for the fermion hard- 
core potential (the potential that gives the exact Laugh- 
lin wavefunction at the lowest Landau level). More re- 
cently, extensive exact diagonalization studies on small 
systems have provided strong evidence that the ground 
state of a partially filled higher Landau level is indeed a 
charge density waveta. The UCDW state is anisotropic 
and it is feasible that its low energy transport properties 
are anisotropic as well. 

There have been suggestions, however, that, for a 
Coulomb interaction, the ground state is a more com- 
plicated structure than a UCDW. A renormalization- 
group analysis by MacDonald and Fisher (MF) sug- 
gests an instability of the UCDW against formation of 
a (highly anisotropic) Wigner crystal in the presence of 
backscatteringu. Crystallization was also found in nu- 
merical Hartree-Fock calculations of Cote and FertigB. 
However, Cote and Fertig found the interactions of the 
modulations among different stripes to be so weak as 
to render the effects of the Wigner-crystal instability 
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immeasurable. The long wavelength behavior oL.the 
anisotropic crystalline state, derived in that workfj, is 
indistinguishable from that of a UCDW state, and, for 
practical purposes, one can ignore backscattering as neg- 
ligibly small and consider the ground state to be a 
UCDW. Thus, throughout this paper we will make the 
assumption that the ground state is well described as a 
UCDW. 

To go beyond the ground state properties one needs 
a way of describing excitations around the ground state 
and effective interactions between them. This is the main 
objective of this paper. The low energy excitations can be 
described as smooth "elastic" deformations of the UCDW 
ground state. We represent these excitations as coherent 
states and derive an effective interacting Luttinger liquid 
model for them, complete with analytical expressions for 
the interaction matrix, on the basis of a first-principles 
Hartree-Fock calculation. 

The interaction matrix of the Luttinger liquid model, 
which we calculate, can be expanded for long wavelengths 
to get what essentially are the low energy elastic param- 
eters of the system. These quantities can be used as 
inputs for calculating a great number of physical quan- 
tities. In Ref. [ll] the elastic parameters were used as 
an input for a calculation of transport properties of the 
system. In Ref. ||, MF use these parameters as inputs 
for the renormalization-group analysis of backscattering. 
The parameters are also used in Ref. |l2] to study the 
Kosterlitz-Thouless like dislocation unbinding transition 
of the striped phase. The calculation of the elastic pa- 
rameters in that reference agrees with our work in the 
long wavelength limitlij. (Our work, however, is also 
valid at higher wavevector). Similarly, in Ref. |l3| these 
parameters are used as inputs for an RG calculation of 
the effects of d isor der on the striped statefij. Finally, 
in Sec. VI and VII of this paper, we will use these pa- 
rameters in calculations of the dynamical properties and 
the I-V spectrum for tunneling of electrons into a striped 
state. We certainly expect that many more applications 
of these results will be found. 

In Sec. H we will outline the expected form of the low 
energy Luttinger liquid theory - including a detailed dis- 
cussion of symmetry considerations. We will also make 
an explicit connection between_the form of the Luttinger 
liquid theory proposed by MFEj and the low energy phe- 
nomenological theory proposed in Ref. ll| . 

Section III will be devoted to the description of the 



find the parameters of the low energy theory of Ref. [n} 



low energy excitations in terms of semi-classical coherent 
state operators. We begin by discussing the Hartree- 
Fock ground striped ground state. Deformations of this 
ground state are described as coherent excitations of el- 
ementary density operators. We will show how our co- 
herent state description naturally leads to a quantized 
Luttinger liquid Hamiltonian. 

In Sec. IV we will calculate the energies of low energy 
coherent state excitations and hence obtain the full inter- 
action matrix of the Luttinger liquid theory. In Sec. [v|, 
we expand the interaction matrix for long wavelengths to 



Sections VI and VII will give examples of applications 
theory and its parameters to the 



of the Luttinger liquic 
dynamics of and tunneling into the striped state 



II. FORM OF THE LOW ENERGY THEORY 

In this paper, we assume that the UCDW ground state 
is a set of straight electron, (or hole) density stripes, con- 
tinuous in the x-directiorJj. Throughout the paper, each 
stripe is labeled with integer indices / or J. The edges 
of the stripes are labeled a or f3 that take the values R 
(or "+1") for the right edge and L (or for the left 

edge. The requirement that the UCDW state minimize 
the energy fixes the distance between the centers of the 
stripes, a. All our microscopic analysis is restricted to 
the partially filled Landau level. The only contribution 
of the lower, completely filled, Landau levels to our model 
is screening of the Coulomb interactionsli-l 

The general form of the Luttinger liquid model 
has beeji written down by MF on phenomenological 
groundsu. In their work, MF assume that the low en- 
ergy excitations are small displacements of the edges of 
the stripes. The dynamics of the edges is represented 
by a quadratic Hamiltonian for the displacement fields, 
u a (q). MF argue that the displacement fields are asso- 
ciated with the chiral currents on stripe edges, p a (q), as 
M a(q) = ot2itl 2 and the Hamiltonian for the cur- 
rcnts-ias a form very similar to that of the displacement 
fields^: 
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(1) 



The interval of integration over q x is cut off by 27r//; 
q y is integrated over the interval [— ir/a; ir/a]; D a ,a(q) 
is the effective interaction matrix. Except for a prefac- 
tor of l/a, this is exactly the Hamiltonian given by MF. 
(The prefactor appears because of the difference between 
our Fourier transform conventions.) This classical the- 
ory is quantized by imposing Kac-Moody comm utation 
relations onto the currents pi a (q). (See Sec. IIIB below.) 

As MF point out, the general form of the elements of 
the interaction matrix, D Ql g(q), .ap. constrained by the 
symmetries of the "striped" stateotL3. Symmetry consid- 
erations become more obvious if we introduce symmetric 
(S) and anti-symmetric (A) modes that for each edge are 
defined as 



PA« = PR.,^ a ^ + P L,qer iq » a! T , 



(2) 
(3) 



where v* = v — [u] is the fractional filling of the par- 
tially filled Landau level. The phases e ±l iy aE r enforce 
the particle-hole symmetry of these modes at v = 1/2. 
These factors are also natural in the sense that the right 
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(left) edge is displaced a distance ±av* /2 from the center 
of the stripe, and thus an additional phase is included in 
its q y Fourier transform. 

The Hamiltonian for the symmetrized excitations be- 
comes 

J2 



H= ^J dq [^(-q)K q | 2 + £>A(-q 



-q)|PA,q| 



ImipA^ps^jlmie^' D LR (-q)} 



(4) 



where the symmetric and antisymmetric coupling con- 
stants relate to the interaction matrix D Q| g(q) as 



^(q) = ^e^^( a -«i)^(q) 



(5) 



This change of variables separates the two components 
of any excitation around the "striped" ground state: the 
displacement of the guiding center and the change in 
width of the stripes. 

The displacement of the guiding center of each stripe 
is described by the symmetric component, ps,q- In this 
case, the excitations on the opposite edges of each stripe 
are in phase. This is equivalent to saying that the width 
of each stripe stays the same, but the guiding center of 
the stripe fluctuates around its equilibrium position. As 
was mentioned by MF and others, the form of the elas- 
ticity for this mode should be like that of smectic liquid 
crystals, £>s(q) oc c\q^ + i.e. there is no q\ con- 

tribution. To see why the q% term should be absent, we 
consider a long wavelength symmetric perturbation in 
the (^-direction. Locally, such an excitation looks like a 
"shear rotation" as shown in Fig. 



the stripes (for fixed amplitude of perturbation) is pro- 
portional to q% (since the angle of the shear rotation with 
respect to the unperturbed state is proportional to q x ). 
However, the initial stripe spacing a is already chosen to 
minimize the ground state energy of the system Eqs so 
that dEas/da = 0. Thus, the leading term of the energy 
cost of the shear rotation is proportional to the square of 
the change in spacing, or as q^.. It also must be noted, 
that the rotational symmetry can be broken by disorder 
or the geometry of the sample and the q% component may 
appear, changing various properties of the system, such 
as RG scaling or tunneling. 

The anti-symmetrized field, pA.q, on the other hand, 
represents the local change in the width of each stripe. 
At distances large compared to the period of the UCDW, 
a, the stripe width fluctuations look like an electron den- 
sity wave, so that the matrix element Da (q) must be pro- 
portional to V(q), the appropriately screened Coulomb 
interaction. The leading corrections to this singular 
Coulomb form of the long wavelength density-density in- 
teraction will be called the inverse compressibility x 1 - 
This inverse compressibility is simply related to the en- 
ergy change for changing the density if the Coulomb piece 
is compensated by a background charge that leaves the 
system locally neutral on scales larger than the stripe 
spacing. 

Using only these kind of symmetry arguments and the 
similarities between the UCDW states of the 2DEC and 
the smectic liquid crystals, Fogler and VinokurEJ (FV) 
proposed an equivalent hydrodynamic Hamiltonian. (A 
similar elastic description is written in Ref. Ex- 
panded for long wavelength excitations to second order 
in the fluctuations, the Hamiltonian of FV may be writ- 
ten as H FV = H + H K , where H K is the kinetic energy 
term, and 
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FIG. 1. A long wavelength symmetric excitation lo- 
cally behaves like a shear rotation. Such a rotation changes 
the distance between the stripes, so that the new distance 
a' = acosO ps ail — 6 2 /2), where 6 is the angle of the local 
rotation. The angle 6 is proportional to q x , so that a — a' oc q 2 . 

Although a uniform rotation costs no energy, this local 
shear rotation also has the effect of changing the perpen- 
dicular spacing of the stripes. The change in spacing of 



H =\j^mi+Kqt)\u^ 

+ (V(q) + X _1 )|5n q | 2 + 2imC'q y u cl 5n^}. (6) 

The hydrodynamic fields are n q , the coarse-grained com- 
ponent of electron density, (mn is the mass density) and 
M q , the displacement of the guiding center of each stripe. 
The field (5n q is the displacement of n q from its equi- 
librium value no = ^jz- The coefficients Y and K are 
the compression and bending elastic moduli, V(q) is the 
appropriately screened Coulomb potential, and x 1S the 
compressibility. Here, the inverse compressibility, as dis- 
cussed above is just the long wavelength density-density 
interaction once the Coulomb interaction is "removed" 
by compensating with a neutralizing background charge. 
Since the ground state energy is always measured with 
such a neutralizing background, the inverse compressibil- 
ity can be written as the change in energy when the stripe 
width is changed by a small fraction while the distance 
between the centers of the stripes is kept constant, 
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with Eqs the energy of the striped ground state. The 
constant C, above, is the measure of the change in the 
local period of the stripes with the change in local elec- 
tron density. The condition of equilibrium derived from 
the Hamiltonian gives a relationO 



C = - 



Y 0(a) 
ma dna 



(8) 



This is obtained by requiring that for a given long wave- 
length variation in n(r), with q x = and q y — ► 0, 
the energy should be minimized with respect to u (i.e. 
5H I Su = 0), when the change in the local lattice constant 
(5a = adyu) has the value required by local equilibrium 
(6a = Snd(a) /dn ). 

Since the field u q in the Hamiltonian is the displace- 
ment of the guiding center of each stripe, it corresponds 
to the symmetrized field ps,q 



u a = 2nr — — . 
q 2 



(9) 



The other field in the Hamiltonian, 5n q , is the deviation 
of the coarse-grained electron density from the equilib- 
rium value, no', its counterpart in the Luttinger liquid 
description is the anti-symmetrized field PA,q- Smearing 
PA,q over the stripe width leads to the correspondence 



8n c 



(10) 



Thus, the first three terms of the hydrodynamic Hamil- 
tonian exactly correspond to the three terms in the Lut- 
tinger liquid Hamiltonian as represented in Eq. (^) . The 
last two terms are introduced in Ref. [Il] for bookkeeping 
purposes only. 

The kinetic energy term Hk disappears when one im- 
poses a strong magnetic field and projects onto a single 
Landau level. The projection, however, has the effect 
of changing the commutation relations, so that the vari- 
ables u and 5n no longer commute with each other. (The 
changed commutation relations were not explicitly stated 
by FV, but were used implicitly in their equations of mo- 
tion). 

FV included an additional counter-term in their Hamil- 
tonian of form 



\[d y i 



(11) 



This counter-term is necessary when terms beyond sec- 
ond order are included in the Hamiltonian, in order to 
assure that rotational symmetry is maintained and that 
(Vk) = in equilibrium. Here, we only need the Hamil- 
tonian to second order, so the counter-term is absent. 

For convenience, we adopt the notation of the hydro- 
dynamic Hamiltonian to write down the expected form of 



the Luttinger liquid interaction matrix D(q) in the long 
wavelength limit 



£>s(q) = 
£>A(q) = 
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In the course of the paper, we will develop a micro- 
scopic Hartree-Fock description of the "striped" state 
and see the Hartree-Fock theory exactly confirm the phe- 
nomenological expectations of the form of the Luttinger 
liquid (and, thus, hydrodynamic) Hamiltonian. The mi- 
croscopic description will also give us quantitative results 
for all the physical parameters of the effective Hamiltoni- 
ans above as well as the necessary commutation relations 
to derive the correct dynamics of the system from the ef- 
fective Hamiltonian (Eq. ||) of FV or the Luttinger liquid 
Hamiltonian (Eq of MF. 

In Section n| below we show how our coherent state 
description of low energy excitations naturally leads to 
a quantized Luttinger liquid Hamiltonian that has the 
form of Eq. (|l|) and can, in principle, be used to obtain 
the anharmonic terms not discussed in this paper. We 
will in fact find that the Luttinger liquid Hamiltonian is 
precisely equivalent to Hartree-Fock within the low en- 
ergy subspace defined by linear combinations of p a (q)\G) 
(\G) is the ground state). 

In Sec. |y| we will calculate the energies of low en- 
ergy coherent-state excitations and, in Sec. obtain the 
coefficients Y, K, C, and x using the Hartree-Fock ap- 
proximation. This is the main quantitative result of the 
paper. A parameter similar to Y has actually been cal- 
culated by MF. 

The parameters Y, K, C, and x fully define the in- 
teraction matrix D a ^(q) in the long wavelength limit, as 
well as the hydrodynamic Hamiltonian (Eq. (^)). These 
parameters have also been calculated in Hartree-Fock by 
Wexler and Dorseyli3, and agree with the long wavelength 
limit of our results. Our results, however, are more gen- 
eral and are valid at smaller wavelengths. Thus, we will 
have presented a step-by-step derivation of the effective 
interacting Luttinger liquid Hamiltonian with all param- 
eters on the basis of a microscopic Hartree-Fock calcula- 
tion. 



III. MICROSCOPIC DESCRIPTION OF 
EXCITATIONS 

A. Hartree-Fock Ground State 

To derive a low energy effective theory from a micro- 
scopic Hartree-Fock theory, one has to find a way of treat- 
ing low energy excitations within Hartree-Fock. But first, 
as a reference point, one has to construct a ground state. 
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A Hartree-Fock state is a single Slater determinant 
made up of one-electron wavefunctions. Since the 
ground state, |G„«), is assumed to be a UCDW, trans- 
lationally invariant in the x-direction, the one-electron 
wavefunctions of choice are the solutions to the free- 
electron Schrodinger equation in the Landau gauge, A = 
(-By, 0,0): 



4>l{v) 
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c <My - k Xn i 2 ) 



where L x is the length of the system in the ^-direction 
(assumed to have periodic boundary conditions), k xn — 
^-n,n = 1,2,.., A, and A is the number of electrons 
needed to fill a Landau level. The pseudo-momentum k x 
controls the position of the wavefunctions (the wavefunc- 
tions are centered around k x l 2 in y-direction) and labels 
the states. A striped state of non-interacting electrons 
with filling fraction v* and distance between the stripes 
a can be created by filling all the wavefunctions with 
k x € [na/l 2 , (n + v*)a/l 2 ] (n is an integer) and leaving 
the others empty. 

For future use, we construct a vector of the com- 
plete set of the one-electron wavefunctions in the Lth 
Landau level: 



Ipk 



(13) 



We restrict the Hilbert space of the one-electron states 
to the highest partially-filled Landau level, L, assum- 
ing that the presence of lower completely filled Landau 
levels manifests itself only in screening of the Coulomb 
interaction!!.! 

We can define an anti-symmetrizing operator A that, 
when acting upon ^l, creates a Slater determinant of 
the set of the one-electron wavefunctions comprising \& £. 
This Slater determinant would represent a Hartree-Fock 
wavefunction for a completely filled (v* = 1) Hartree- 
Fock state. 

To represent the striped ground state, \G V *), however, 
we need to make a Slater determinant of a subset of 
the one-electron wavefunctions that correspond to the 
filled states. By analogy with A, we can define an oper- 
ator A a ,v* that creates such a Slater determinant for the 
striped state with filling factor v* and distance between 
the stripes a. When it acts upon ^ l, the operator A a ,v* 
effectively picks out the wave-functions of the filled one- 
electron states and anti-symmetrizes them. We use the 
resulting Slater determinant, A a .v*^L, as the Hartree- 
Fock trial ground state of the striped quantum Hall sys- 
tem. In the spirit of variational Hartree-Fock, the width 



of the stripes, a, is the parameter to be adjusted by min- 
imizing the Hartree-Fock energy of the state, as will be 
shown in Sec. 



The Slater determinant A a ,v'^L with 
the optimal distance between the stripes is our Hartree- 
Fock ground state. 



B. Coherent States 

Having constructed the striped ground state, we can 
proceed to describe low energy excitations around it. In 
this section, we will introduce a method for creating any 
low energy excitation around the striped ground state in 
a way that connects Luttinger liquids with Hartree-Fock 
states. 

We start from the Luttinger liquid description of ex- 
citations. As was mentioned in the Introduction, we as- 
sume that the ground-state UCDW has a striped struc- 
ture in momentum space, with the filling factors in the 
stripes alternating between integers \v] and [v] + 1 . Effec- 
tively though, the filling factors in the stripes are and 
1, since we restrict all microscopic analysis to the top 
Landau level. Thus, since backscattering is ignored, we 
can consider the edges of the stripes to behave like the 
edges of large integer quantum Hall states, with filling 
factor 1. In such states, the excitations around an edge 
of the system are often described by a density operator: 



(><l,. 



k x +q. 



where Cfc x annihilates an electron in a one-electron state 
labelled by k x . The operator p 9x describes a current 
running in one direction along the edge of the sys- 
tem. It obeys the Kac-Moody commutation relation 
{\pqxiPq*'] = QxS(q x + q x ')) that determines its algebraic 
behavior and dynamics. Thus, on one hand, this op- 
erator is defined microscopically, through creation and 
annihilation operators, and, therefore, can be used in the 
microscopic Hartree-Fock analysis. On the other hand, it 
can also be described through its algebra as part of a Lut- 
tinger liquid description of the system, thus providing a 
connection between Luttinger liquids and Hartree-Fock. 

However, we cannot use this operator to describe ex- 
citations localized to a particular edge of the stripe since 
p qj . acts across the entire system. When applied to the 
striped ground state, it would create the same excitations 
on all edges of all stripes. To treat the excitations on a 
particular edge of one of the stripes, we need to local- 
ize the operator to that edge. We do so by introducing 
an envelope function or "wavepacket" function wi a (k x ). 
The function wi a (k x ) can have an arbitrary form as long 
as it equals unity in the vicinity of the edge al and the 
different wi a functions have disjoint support — meaning 
only one wi a (k x ) can be nonzero for any given k x (and w 
should not have any unintegrable singularities). It is im- 
portant to note that all the physical quantities calculated 
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in this paper are independent of the shape of wi a (k x ) 
away from the edges, and one never actually needs to 
worry about the specific form of wi a as a function of k x . 

With this "w-wavepacket" , we define a new density 
operator that affects only the given edge al: 

Pal,q* = J2 Wla ( kx + fe/ 2 ) c L+<fc C fcx> ( 14 ) 

where w a j(k x ) has to be unity within a range of ±q x 
around the edge al to ensure unitarity of the operator 
U e aI q ^ defined below. 

The new density operator obeys the Kac-Moody com- 
mutation relation 

[PaJ.fe , P0J,qA = aq^udapS^ + q x '), (15) 

just like the ordinary density operator p qx , but localized 
around the given edge. Thus, it is this operator p a j „ x 
that enters the Luttinger liquid Hamiltonian (Eq. (jlj)) 
but is also defined in terms of microscopic electron cre- 
ation operators (Eq. Q)), thus allowing us to make a 
connection between microscopic Hartree-Fock and the 
Luttinger liquid Hamiltonian. 

To complete our description of low energy excitations 
around the striped ground state, we construct an opera- 
tor 

U e aIqx = e^'-'.fa+^^.faW (16) 

that creates a state 

\aI,q x ;e) = U e al! JG) (17) 

when it acts upon the ground state. (\G) = \G U ») fur- 
ther on to simplify notation.) The operator U^j is con- 
structed by analogy with coherent states of a harmonic 
oscillator. The operator p a i,q x acts as a creation oper- 
ator and e as the amplitude of the localized excitation 
that p a i,q x creates. 

One can check that e does in fact play the role of the 
amplitude of the localized excitation by calculating the 
classical expectation value of the electron density opera- 
tor, pj,p,q^i, in this state, 

{al,q x ;e\ppj^\al,q x ;e) = 
aSijS al3 (ieq x S(q x - q x ) - ie*q x S(q x + q x ')), 

for a small |e|. This relation implies that the coherent 
state is indeed a state that corresponds to a classical 
electron-density excitation on the given edge al. This 
density excitation is real and sine- like (cosine-like), with 
wavevector of magnitude \q x \ if e is real (imaginary); the 
amplitude of the excitation is proportional to |e|. 

Since we substitute what would be a state, in which 
a complex chiral current runs along the edge al, with 
two states with real (cosine and sine) density excitations 
around that edge, we restrict q x to q x > to avoid over- 
counting. 



One may notice that excitations with q x — cannot be 
created in this manner. However, this is not a limitation. 
In the absence of backscattering the number of electrons 
(holes) in each electron (hole) stripe is conserved. It is 
easy to see that excitations with q x — violate this con- 
straint and thus cannot be present in our approximation. 

Making localized excitations around the edges of dif- 
ferent stripes, we can create any possible low energy ex- 
citation around the striped ground state. We describe a 
set of excitations along the edges al in terms of their x- 
wavevectors q x ,ai and their complex amplitudes e a i. The 
most general form of the coherent state \{q x , a i] £ai}) is 
one in which there are different density excitations on 
every edge. Since the operators U e aI q ^ commute for dif- 
ferent al we can create excitations of multiple edges just 
by successively applying these operators - one for each 
deformed edge. 

In this work, however, we will only be concerned with 
states in which all edges have the same x-wavevector q x 
(since it can be easily shown that excitations with differ- 
ent q x do not interact with each other to quadratic order 
in the deformations). We thus define an operator 

U(q x ; {£*/}) = II^X ( 18 ) 

la 

which, when applied to the ground state, yields the state 

\q;{e aI }) =U(q x ;{e aI })\G); (20) 

which has amplitude |e a /9a;| on stripe al and all edges 
have the same q x . 

We will construct and calculate the energy of small- 
amplitude coherent states in both the Luttinger liquid 
language and in Hartree-Fock language. Equating these 
energies will allow us to determine the full form of D a p{q) 
of the (previously) phenomenological Luttinger liquid 
Hamiltonian. Once such a Hamiltonian is written down, 
it is equivalent to Hartree-Fock in the sense that any of its 
matrix elements within the subspace defined by p a i iqn . \ G) 
are precisely equal to the corresponding matrix elements 
calculated in Hartree-Fock. 



C. Energy of Coherent States in Luttinger Liquid 
Language 

The localized density operators p a i.q x obey Kac- 
Moody commutation relations and describe currents run- 
ning along given edges al. These fields, therefore, are 
exactly the fields we need to construct the effective Lut- 
tinger liquid Hamiltonian (Eq. (fflh) 

H = ttI 2 / d 1xa(3p a i^D af3 (~q x ,I - J)ppj,-q x 

a,p i, J 

(21) 



G 



Here the Hamiltonian is Fourier transformed in the re- 
direction, along the stripes. 

Semiclassically, the energy of a given state equals the 
expectation value of the Hamiltonian in this state. The 
excitations of interest to us are the low energy excita- 
tions that, in the previous subsection, we represented by 
the coherent states and whose dynamics the Luttinger 
liquid Hamiltonian describes. Given a coherent state 



\%c\ {Sal}), the classical expectation value of this energy 
in this state is then 

E= (q x ;{e aI }\H\q x ;{e aI }) (22) 

If the interaction matrix, D, is known, for small values 
of |e a i|) we can expand the energy to second order in 
|e a i| to obtain 



E = E + irl 2 E / dq' xl 5D lS (-q' xl M - TV) X 

aI,0J jM,SN 

[(G\(E a lp a i tqx +e* a i{p a I,q x f)p 1 M,q' x P5N-q' x { e PJP0J,<l* + e pAPPJ,q*f)\ G ) 
-^{G\{e a lp a l, qtc + e* aI {p aI ^){epjppj^ + £* j{Pl3J,qJ { )p 1 M,q' x P&N,-q' x \G) 

-^( G \PiM^PSN- q ^{e aI p aI ^ +e*/(p Q /, g J t )(e/3jP/3j, 9x + s* 0J {ppj,q x ) ] )\G)}, (23) 



where the constant Eq = (G\H\G) can be considered to 
be the ground state energy. This ground state energy is 
a large additive constant that is the same for all excita- 
tions. It appears because the coherent state describes the 
entire excited "striped" state and not just the Luttinger 
liquids on the edges of the stripes (the set of Luttinger 
liquids can be thought of as the excited "striped" state 
minus the "striped" ground state). Given values of the 
coefficients ei a , Eq. (|2^) can be evaluated using the com- 
mutation relation of the density operators p a i q given by 
Eq. @. 

We now create an excitation of wavevector q = (q x , q y ) 
by taking e a i to be given by 



e Ia = rje 



(24) 



where rj is the overall amplitude of the excitation. We 
recall that the complex phase of the parameters e de- 
termine the phase of the excitation along the edge - 
representing a sine excitation if e is real or a cosine if e 
is imaginary. Thus, when the constant 6s = we have 
a symmetric (smectic) excitation; when 9s — tt we have 
an antisymmetric excitation; and for #5 = tt/2 we have 
a "mixed" excitation, in which there are a cosine excita- 
tions on one side of all stripes and sine excitations on the 
opposite side. Inserting this form into Eq. (^3|) yields 

E = E + ^^x 
2 a 

D LL (q) + D RR (q) + 2Re{e t9s+ ^' D RL (q)} 

yielding the following three energies for symmetric, anti- 
symmetric, and "mixed" excitations 

AEs = fofeQ 2 L X N ^ eift ,4 (a - ff ^ (q) 



AEa 



2 a 



(fjqj) 2 L X N 
2 a 



a/3 



5^a J 8e i *' Tr(«-/ J )D a/J (q) (25) 



AE 



M 



(■qqj) 2 L X N 



[Dll(cl) + D RR (<i) 



2Im{e tq * c 



If the interaction matrix D is known, we can obtain 
the energies of any small excitation using the Luttinger 
liquid description. Conversely, it is clear that, if the en- 
ergies E are known for arbitrary (small) e, the matrix 
D may be obtained. We shall calculate these energies in 
the Hartree-Fock approximation for the coherent states 
that we can describe microscopically as well as through 
the density operators (Sec. IV). 



IV. ENERGY OF COHERENT STATES 



A. Coherent States in Hartree-Fock Language 



To find the energies Es, Ea and Em (Eq. fl25|)) from 
first principles, we need to describe the low energy excited 
states within Hartree-Fock. This also can be done using 
the coherent states, since the formalism can be easily 
translated into the Hartree-Fock language. To determine 
the energies, it will suffice to consider the states where 
q x is the same for every edge. Within the Hartree-Fock 
approximation, the operator U(q x ; {s a i}) is a rotation in 
the Hubert space of the wavefunctions {V'L,fc x „} of the 
form 



U{q x ;{e a i}) 



(26) 



a/3 



where R(q x ; {e a j}) is a matrix whose elements are given 
by 

Rk x ,k' x (,Qx', {£«/}) = 

= S K ,k x -q x W(k x - y ) + S K , kx+qx W*(k x + |), (27) 
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where 



W{k x 



E< 

a/ 



xlWalikx) 



(28) 



and w a i{k x ) is the wavepacket envelope function defining 
the operator pi a in Eq. (|l4|). 

The rotation U(q x ; {e a i}) is a unitary operator that 
preserves orthogonality and completeness of {?/>l,/c x „}, as 
long as the w-wavepacket is unity over a range of ±q x 
around the edge. The Slater determinant A a ,v* ^ q "' , in 
which iffq^ 1 — U(q x ; {e a i})^L and all the states are filled 
in a striped fashion, represents the excited coherent state 
\lx't { £ ai}) in the Hartree-Fock approximation (Fig. ||). 




FIG. 2. Top figure is a schematic representation of the 
U{q x ; {e a i}) transformation. Each single-electron wavefunc- 
tion, represented by a solid line here, is transformed by 
U(q x ; {£«/}). Then a Slater determinant, Aa,v* ^q" 1 , is cre- 
ated, in which the new single-electron wavefunctions are filled 
in a striped fashion, just like in the unperturbed case. It 
is important to note here that, using the transformation 
U(qx\ {Ed}), we can also create states that have different exci- 
tations on different edges of the system. The electron-density 
profile in the bottom figure illustrates this point. Here, an 
excitation is created on the left edge of one stripe (we used 
the L=0 single-electron functions for simplicity). 

Later in this Section, we will calculate the energy of 
the the Slater determinant A a ,v*^ e q " 1 for e a i correspond- 
ing to symmetric, antisymmetric, and mixed excitations 
(See Eq. ^i]). In this way we will find the energies AE$, 
AEa , and AEm and read off the full interaction matrix 
of the Luttingcr liquid theory, D(q). 



V mnst = / ^ m (r)VL f ,(r)Vb(r-r')Vfc.. (0^1,(0 



2tt 



27T 



P), (29) 



where 



V eS (q x ,q y ) = V L {q) (l L (jf)) C '~ 



(30) 



q = yq x + qy and L^{q) is the Lth Laguerre polynomial. 

The screened Coulomb potential, Vz,(r), is related to the 
bare Coulomb potential, V^(r) = 2ne 2 /(nq) with k the 
background dielectric constant as 



V L (q) 



2ne 2 1 



q Ke(q) ' 



with the correction ter: 



e(q) = l+ [1- J$(qR c )], 

qaB 



(31) 



(32) 



arising from screening of Coulomb interaction between 
electrons in the partially filled upper Landau level by 
the electrons in the completely filled lower Landau lev- 
els. Here as — h 2 n/m*e 2 is the Bohr radius and m* is 
the effective electron mass in the system. Throughout the 
paper, J M (x) is the Bessel function of order /j,. The form 
of the dielectric constant in Eq. (|3^) is obtained through 
a semiclassical approximation and is valid for L » 1 . A 
fuller expression for the effective dielectric constant that 
is valid for small L is given in Ref. [l4|. 

The Hartree-Fock energy of the ground state is there- 
fore 



Egs = 



1 L 2 



in in n >i * in ii n in 



0, 



(33) 



where k aI = a/ 1 2 (I + av*). The first term of Eq. ( |33j ) is 
the direct component of the ground state energy, Eq S , 
and the second term is the exchange component, Eq S . 
The superscripts D and X will label direct and exchange 
energies henceforth. After substitution of Eq. (|2^) into 
Eq. (|3^) r ajid some mathematical manipulations, Eq. (|3J 
becomealj 



B. Ground-State Energy 

In the space of the one-electron wave- functions {ipk xn }, 
the Coulomb potential screened by the lower, completely 
filled Landau levels, Vi(r),til takes on a matrix form 



1 L x Na ^ 2%n. ( sin(7w*n) x ' 
^— ' a \ nn 



U{q) 



2 2ttZ 2 
2tt/ 2 



H^Xff(^p) ; 



(34) 
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where the ratio L x Na/2irl 2 is the number of states in the 
partially filled Landau level and V^g(0, 0) is assumed to 
be 0, to account for the uniform positive background. 

The width of the stripes is fixed by minimization of the 
ground state energy (Eq. (|34|)), i.e. by the requirement 
that d a E GS = 0. 



C. Hartree-Fock Energy of Excitations 



Following the procedure given in Sec. Ill, we represent 
the excitations by coherent states. Within the coherent 
state formalism, any low energy excited state turns out 
to be conveniently represented by a single Slater deter- 
minant, A a , u *& q f, where = U(q x ; {e aI })f (the 
Hilbcrt-space vector ^ is defined by Eq. (|l3|)), that can 
be used as the Hartree-Fock wavefunction. The trans- 
formed single-electron states arc 



i>k x ( 

+ iW(k x 
- W{k x - 
-W*(k x 



l-\W(k x + ^)\ 2 -\W(k x ~^)\ 2 ) (35) 



2 
qx 
2 



Hk^+q* +iW*(k x 



2 

q x 



)^k 



)W(k x + 
)W*{k x 



3g a 



)ipk 



to second order in the small amplitude (W or e a j). To 
the same order, the matrix element of the renormalized 



Coulomb potential in this set of states is 



V„. 



V L (r - r')4, m (r)r kxn (r)^, (r')^ t (O 

*+£[ 



b=±l 



+ W(k xm + b^-)W*(k xn + bY)V(m+bq m )(n+bq x )st + 

q x q x 
+ W(k xs + b—)W*(k xt + b—)V mn ( s+bqm )( t+bq;c) + 



+ W(k xm + b^)W*(k xt 



b^Wim+bq^nis-bq^t + 



W(k xn - b^-)W*(k xt + 6y)F m(n _ fc(?a;)s(t+b(?x) ] 



It is easy to see that, because of the (^-function in Eq. (EQ), 
the last two terms of the expansion of ipk^ (Eq. (pq)) 
do not contribute to the matrix element V mn st until the 
fourth order in W. Effectively, we need to expand tp^ 
only to first order in W and normalize it. 

The Hartree-Fock energy of given excited state is calcu- 
lated the same way as the ground state energy (Eq. (j||)) 
but with (V mmnn - V mnnm ) substituted for {V mmnn - 
Vmnnm)- Slightly tedious but straightforward computa- 
tion leads to the Hartree-Fock energies for the symmetric 
(AEs), antisymmetric (AEa), and mixed (AEm) exci- 
tations above ground state energy. If we define two func- 
tions 



F(k x , ky) — 



V cS (k X ,ky) 



2irl 2 
4gin 2 ( .fc 2 ^fc JL ' 

fay) 



dpdu 



cos(pl 2 k x ) cos(ul 2 ky)V e g(u,p), 



we get a greatly simplified expression for the energies 



AE S 
AE A 
AE Mo 



2M 



,L X N 
2na 



(gxQ 2 ^ sin 2 (imv* ) ( A s (q x , 



2\ v 



L^Na 



2 "as 



27m 



/ 7\2V^/ 9/ *\ a / 27T7J Z7T71 

\1xly 2^( cos t 7 ™* 7 ) A s{qx, Vq y )F{q x , 



, „, 27m . . , 27m, „. 27m,, 
l )F{q x , + q y )- A s (q x , )F(0, )), 



2nn 



= \V 



2na 
L X N 



27m, 



(36) 
(37) 

(38) 



27m, 



1y) 



AE 



2na 
AE S + AE A 



(q x l) 2 ^ Bm(2imv*)A a (q X) 



2im 2~Kn 
\-Qy)F{<lx, \-q y ), 



sm\^nu*)A s {q Xl )F(0, )), (39) 

a a 



(40) 



AE 



Here we used exactly the same coherent states to represent the three basic excitations, as we did in Eq. ( p5|) to obtain 
the general form of AEs, AEa, and AEm from the Luttinger liquid Hamiltonian (note that it is AEm and not 
AE M that is proportional to Im{e i9 « cuy * D RL (q)}Y Using Eqs. © and (|), we can now simply read off the elements 
of the interaction matrix D(q) from Eqs. (Bq)-(ptJ): 



n 



■ 2/ * \ a i 27m 27m 

sm [jxnv )A s (q x , h q y )F(q x , 

a a 



\ ■ ii *\ a / 27m, 27m, 

q y ) - sm [irnv )A S (q x , )F 0, 

a a 



o, * \ a / 27m 27m 2 , , 2nn 2im 
cos (imv )A s (q x . h q y )F(q x , \- q y ) ~ sin (irnv )A s {q x , )F(0, ) 

CL CL G/ CL 



, (41) 
, (42) 
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l m{e W D RL ( q )} = -—y / sin(2nnu*)A s (q x , + q y )F(q x , +q y ), 



(43) 



These expressions for the Luttinger liquid interaction 
matrix D(q) are valid for any wavevector q y and for small 

Now, we can expand Eqs. (^)-(^) for q x l and 
q y l << 1 and re- write them in terms of the hydrody- 
namic parameters Y, K, \i an d C, thus completing 
the connection between Luttinger liquids and hydrody- 
namics and obtaining all the hydrodynamic parameters. 



V. PARAMETERS OF THE HYDRODYNAMIC 
THEORY 

Expansion of Eq. ( |38| ) gives the expected form for the 
symmetric component of the effective interaction 



Ds(q)="(^+^)> 
where the coefficients Y and K are 



(44) 



Y 



K = 



— 77 sin 2 (7w*7i) 

7T0 2 ^ 

n 

1 ^ sin 2 (7r^*n) 
16tt ^ (7m) 2 

(27m) 3 Z 4 



^^ eff( o, g) U + ^ 



d ^r ^^^(th 2 ™^ 27m 2 



2n c V a a 



M 9 ^eff (0, *) U - ^ W ff (0, ff ) | ^ 1 + 3 f — ) 



3a 3 
(27m)? 4 



4a 



2tt 



v eS ^m 4 [M(— i 2 ) + Mi— i 2 ) 



(45) 



(46) 



Note there there is no g 2 term in expansion of the 
energy of the symmetric excitation as follows from sym- 
metry arguments for a smectic excitation (Fig. [I]). For 
interactions with circular symmetry, the g 2 -term is pro- 
portional to d a Ecs that vanishes when we minimize the 
ground state energy with respect to the stripe-width, a. 

The effective antisymmetric component of the interac- 
tion (from Eq. ([39])) behaves as 



£>A(q) 



1 V(q)+X 
2-wl 2 



2tt 



(47) 



in accord with the phenomenological arguments given in 
Sec. ||. The constant x _1 turns out to be 



= cos(2irnv*) 



(48) 



-2-kV 



|^eff((U)£Jo(^ 2 ) 



2-kV 



as expected from Eqs. (@) and (|J). If the 2DEG 
is screened by a metal plate placed a distance d away 
from the surface of the gas, then V(q) in Eq. ( (47|) is re- 
placed by V(0) = lim q ^oV(q). The screened Coulomb 



: (1 



hence 



potential has the form V(q) = 
V(0) = 2Tre 2 d/n. 

The off-diagonal term Im{e* ? '°"* D R l(c[)}, given by 
the expansion of Eq. (f40|), is linear in q y when the 
wavevector is small: 



Jm{e^ av D RL (q)} = ^Cq y , 



where the constant C is 
C = — V" sin(27r? 

+1 



2tt/ 3 



d q V eS (0,q) 



(49) 



df . , 27rn 9 , 



(50) 



Note that the constant C and the compression elas- 
tic modulus Y satisfy the relation C = — — --r 2 -, as, 

J ma otiq ' ' 

again, was expected from phenomenological considera- 
tions above (See Eq. (j|)). 



VI. DYNAMICS 

As the first, simplest application of the model, we can 
consider the dynamical properties of the striped state. 

Using the commutation relations of the density oper- 
ators, Eq. (jl|), we can find the Heisenberg equations of 
motion for the density operators 



idtPaq = [p a q,H] 

a T 



(51) 
(52) 



Fourier-transformation and diagonalization of the equa- 
tion of motion gives us the dispersion relation (keeping 
in mind that, by symmetry, D RR (q) = D LL (q)): 
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±- 



2nl 2 



-'I 



yjD RR (q)D LL (q) - D RL {q)D LR {q). (53) 



The expression under the square root is the determinant 
of the interaction matrix, detD(q). Since the determi- 
nant of a matrix does not change with the basis in which 
the matrix is defined, for any basis of density operators 
that we choose for our Hamiltonian, the dispersion rela- 
tion is 



UJ 



, 2ttZ 2 n — . 
± 9xV detD(q) 



(54) 



For example, it is convenient for us to use p A q and psq 
for further calculations. In this basis, the dispersion is 



G e (r) = (V aI (x, 0)*^(x, r)) cx e~ c ^ , (59) 

where 

C(t) = -(((f> aI (x,0)4> a i(x,T)) - ((f> a i(x,0)<j) a j(x,0))). 

(60) 

To get the time-dependent correlations, we follow MF 
and construct an imaginary-time action from the Lut- 
tinger liquid Hamiltonian (Eq. (|l|))u. In terms of the 
bosonic fields <p a i, the imaginary-time action takes the 
form 



Sn = - 



1 f dqdui 



±- 



4irl 



1 



2 J (2tt) 3 



4> a {q, u))M af3 (-q, -w)(j) {-q, -u), 



q x \j D s {q)D A {q) - -(Im{e*W D RL (q)}) 2 , 
a V 4 

(55) 

and we can use the parameters and equations found in 
the previous section to study the behavior of the disper- 
sion more quantitatively. 

In Figure |3|, we plot the positive-^ branch of the dis- 
persion as a function of (q x ,Qy)-, for small values of q x 
and for q y within the Brillouin zone (0, 27r/a); the other 
parameters of the system are specified in the caption. . 
Note that, even though the functions Ds(q), D A (q), and 
lm{e lqyal ' D R L(q)} are not periodic in q y , it is easy to 
show that the determinant, detD(q), is periodic, in any 
basis of choice. "y 

The modes corresponding to the positive and negative^ — 
branches of the dispersion are p+, q and p-, q , respectively, w 
They can be found to be 



(61) 



where 



M a p(q,to) = a5 af3 - h — —D af3 {c\). (62) 

Zira a z 

The minor differences between Eq. (7) of Ref. || and 
Eq. (p^) here are due to the differences between our 
Fourier-transform conventions. 



0.01 



1 

N 



{D A (q)p S q 




(-Im{e^ D RL (q)} - ^/detD(q))p Aq ) (56) 



0.05 0.1 0.15 0.2 

q x /(27C/l) 



0.2 0.4 0.6 0.8 1 



q/(27U/a) 



N 



"In'' 



N 2 



"r>fli(q)} + VdctD(q))^ q ) 
D A (q)(D A (q) + D s (q)) 



(-Im{e 



(57) 
(58) 



VII. TUNNELING PROPERTIES 

The Luttinger liquid model also gives some interest- 
ing results about the properties of perpendicular tun- 
neling into the striped state. Following the usual treat- 
ment of electrons in one dimension, we represent the 
operator for an electron "running along" the edge al, 
^ai, through bosonic fields <p a i, so that ~ e ^ Q/ . 
The bosonic fields are related to the density operators as 
Pai{x) = otd x (j) a i(x)/2Tt. These identifications allow us 
to calculate the electron Green's function 



FIG. 3. Dispersion relation of a striped state with 
v — 4.45 (L = 2, v* = 0.45, i.e. a filling fraction around 
9/2). Coulomb interactions here are assumed to be screened 
only by the lower, fully filled, Landau levels. This screening 
depends on the density of electrons in the system, which is 
taken to be n s = 2.67 x 10 11 cm -2 . In the left plot, the disper- 
sion relation is given as a function of q x for fixed values of q y : 
q y — is the solid line, 0.1(27r/a) is the dotted line, 0.2(27r/a) 
is the dashed line, 0.3(2-7r/a) is the long-dashed line, 0A(2n/a) 
is the dot-dashed line, and 0.5(27r/a) is the grey line. Here, 
the distance between the centers of the stripes that minimizes 
the energy of the striped ground state is a — 6.54L In the 
right plot, the dispersion relation is given as a function of q y 
for fixed values of q x : 0.04(27r/Z) is the solid line, 0.08(27r/Z)is 
the dotted line, 0.12(2tt//) is the dashed line, 0.16(2tt/Z) is 
the long-dashed line, and 0.2(27r/Z) is the dot-dashed line. 
Note that for q x 7^ 0, the curves do not reach u — for 
q y = and 2k/ a but come down to a very small value 
lu — q^l 2 yj K (V (qx) + X -1 ) (that corresponds to the solid 
line in the left plot). 
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The bosonic Green's function is just the inverse of the 
kernel M Q| g(q, uu): 



6 a (q,w)0 /3 (-q, -w)> = M ^q,^) 



(?*') 4 



^detD(q) 



(63) 



Fourier transformation of Eq. ( p3\ j gives an expression for 
the real-space bosonic Green's function that, after inte- 
gration over uj takes the form 



xl{x,Q)<t>al{x,T)) = 

47T J q x 
a f dq D aa {q) a _ M 22dl^p T 



2 J \q x \ ^/detD(q) 



(64) 
(65) 



The first integral vanishes if one assumes that it should 
be taken as a principal value. The second integral can be 
evaluated approximately for large values of r, at which 
our low-energy theory is valid, to give us the low en- 
ergy contribution to the bosonic Green's function and, 
thus, to the electron Green's function, G e (r). It is im- 
portant to note, however, that we should keep to ener- 
gies large enough (times low enough) to be able to ignore 
possible backscattering effects. The energies that we con- 
sider must be higher than the backscattering energy, 1 
where tt, is the backscattering time. 

Note that if the excitations along the edges of the 
stripes were independent of each other, detD(q) would 
equal D 2 Ya and we would get exactly the logarithmic form 
expected of C(t) for a system of non-interacting Luttinger 
liquids. Specifically, we would find G e (t) oc 1/t, as re- 
quired from a set of non-interacting v = 1 edges. How- 
ever, the excitations on the edges of the striped state do 
interact with each other, and the tunneling exponent be- 
haves differently from a noninteracting Luttinger liquid. 

For unscreened Coulomb interactions, the contribution 
of the low energy excitations to C(r) (Eq. (p5|)) can be 
obtained using the approximate low wavevector form of 
D(q) given by Eqs. 



4J, (47), and (Ugh 



C(r) 



2nl< 



Vcl 
Y 



dq 



\qxq v \y/q 



(1 



(66) 



where Vc = 2it1 2 /(Ik) is the characteristic Coulomb en- 



ergy of electrons sc 
Y is given in Eq 



mrated by the magnetic length I, and 
To get the approximate form of 
the interaction matrix, we assumed that q is low enough 
that V(q) » x^ 1 and V(q)Y » (mC) 2 /4 (the expres- 
sion for C is given in Eq. (^p|)). We also ignored the 
q x dependence of the symmetric term of the interaction 
matrix, since it gives rise only to a sub-leading correc- 
tion. The integral in Eq. (66) can be evaluated exactly 
by changing variables to polar coordinates to get 



C(t) 



8T( 
3 



r(i) 



-2 s ) v 2 



■ { YP T) 



(67) 



The prefactor -|r(-|)/r(|) w 8. The time scale 
t = (Yl 2 )/V 2 is very small. For example, evaluating 
Y using Eq. @ for L = 2, v* = 0.4, we find that r 
approximately equals 0.006/Vc and is still smaller for 
higher Landau levels (for L = 3, tq = 0.0038/Vc; for 
L = 6, r = 0.003/Vc, etc.). 

Similarly, for screened Coulomb interactions, we have 



C(r) = 



Til 2 



v(o) + x- 



Y 



\Qxq y \ 



), 

(68) 



where O = ^\/Y{V(0) + x~ r )- I n principle, the in- 
tegral in Eq. ( pq ) can be expressed exactly in terms 
of a generalized hypergeometric function, C(t) = 

-ir 2 l/aQ T), how- 



2 

Til' 2 



v(Q)+x- 



Y 



5 F 3 ({1,1,1},{2,2,2}, 



ever, an approximate but intuitive form may be more 
useful. To analyze the leading behavior of the inte- 
gral as a function of r, we split the integration do- 
main into two: one in which \q x q y \l 2 < l/(^o' r ) (where 
Oo = l/2-\/y(y(0) + x -1 )) an d the exponential function 
in the integrand can be expanded into a Taylor series, 
and another, in which \q x q v \l 2 > l/(f2or) and we assume 
that the exponential is much smaller than 1 and can be 
ignored. The leading behavior comes from the bound- 
aries of the domain, in which \q x q y \l 2 > 1/(S1ot), an d we 
obtain 



7Tf z 



V(0) 



±2Cilog 2 (^^„r). 
Y a 



(69) 



Here the time-scale l/(7r 2 ^f2 ) is proportional to 1/Vc 
with a proportionality constant of order unity. 

It is interesting to compare these results to the tunnel- 
ing exponents of a fractional quantum Hall state at half- 
filling. Unlike the striped state, the fractional quantum 
Hall state is isotropic, but its upper (or only) Landau 
level is also partially filled. As one would expect, tun- 
neling into the isotropic statejis more suppressed than 
that into the striped stateOHB C(r) behaves as r 1 / 2 for 
unscreened Coulomb interactions and r 1 / 3 for screened 
Coulomb interactions. The results are summarized in 
Table I 

The low frequency behavior of the spectral function, 
which can be obtained by taking the inverse Laplace 
transform of G e (r), gives the, behavior of the current- 
voltage relations (I-V curves) .00 The results for I(V) 
are summarized in Table [n|. 



VIII. FURTHER APPLICATIONS AND 
CONCLUSIONS 

The Luttinger liquid model and its parameters derived 
in this paper can be applied to other calculations and 
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studies of various properties of the striped quantum Hall 
state. Thus, for example, Wexler and Dorsey use the 
Luttinger liquid parameters in the long wavelength limit 
to study disclination unbindingE3. They use the param- 
eters as a starting point in their RG analysis and find 
the Kosterlitz-Thouless temperature for the transition 
between orientationally ordered and isotropic phases of 
the striped state. 

MF use the Luttinger liquid interaction matrix at q x — 



to study the st£ 
to backscattering 



jility of the striped state with respect 
. They add a backscattering term to 
the Luttinger liquid Lagrangian and analyze the flow of 
the backscattering amplitudes when u and q x are rescaled 
(while q y serves as a fixed parameter). Their RG analy- 
sis shows that backscattering is always relevant, and, at 
very low temperatures (<10mK), the striped state is un- 
stable with respect to formation of an anisotropic Wigncr 
crystal. We have reproduced their calculation and found 
that, despite some definitional differences, our result for 
the scaling dimension of the backscattering term is iden- 
tical to that of MF. 

We also considered the scaling of the backscattering 
term in the presence and absence of rotational symme- 
try, since in a real sample there are always preferred di- 
rections because of the geometry of the sample and the 
pinning of the stripes by impurities. The assumption of 
rotational symmetry manifests itself in the form of the 
smectic elasticity Ds(q) cx Yq y + Kq x . In the absence of 
rotational symmetry the quadratic term q\ in the elastic- 
ity does not vanish and, in the long wavelength limit, the 
elasticity takes the form: Z?s(q) oc Yq y + K 2 q x . The RG 
analysis can depend on the presence or absence of rota- 
tional symmetry if one chooses to rescale both q x and q y 
(as opposed to treating fixed parameter). How- 

ever, for short-range interactions, the scaling dimension 
of the backscattering term is independent of rotational 
symmetry, i.e. the absence of rotational symmetry does 
not increase the stability of the striped state. It is im- 
portant to note, however, that the rotational symmetry 
might make a difference in stability of the striped state 
with long-range interactions. 

A much more detailed RG analysis of the effects of dis- 
order is given in Ref. [R^. In their paper, Scheidl and von 
Oppen use an elastic theory of the striped states, which 
is very similar to the Luttinger liquid theory, to map 
out a phase-diagram of different disorder and length-scale 
regimesEj. 

The hydrodynamic model, for which we found all the 
parameters by establishing a one-to-one correspondence 
with the Luttinger liquid theory, has been introduced and 
used by FV in Ref. |ll to find dynamical scaling exponents 
of the striped state, 

In addition to these applications and the applications 
described in the previous two sections, the model and the 
parameters will probably find a wide range of other uses 
in the analysis of the striped quantum Hall states. 

In summary, we have derived an effective interacting 



Luttinger liquid model for the low energy excitations of 
the striped quantum Hall state on the basis of a first- 
principles Hartree-Fock calculation. Using our coherent 
state formulation we were able to describe all possible low 
energy excitations (with low q x and any q y ) of the striped 
state, and thus derived the full interaction matrix for the 
Luttinger liquid model, complete with all parameters. 

In the limit of long wavelengths, we have established a 
one-to-one correspondence between the Luttinger liquid 
model and the hydrodynamic model of Ref. [ll]. Thus, we 
have also obtained expressions for the parameters of the 
hydrodynamic model from Hartree-Fock. 

We have also given some examples of applications of 
the Luttinger liquid model and its parameters. In Sec. VI 
we have found the dispersion rela tion for the excitations 
of the striped state and in Sec. VII we have discussed 



the behavior of perpendicular tunneling into the striped 
state and obtained the corresponding I-V relations. Fur- 
ther applications of our results are almost certain to be 
found. 
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u* « 1/2, L > 1 
striped 


v = 1/2, L = 
isotropic 


unscreened 


r l/3 


r l/2 


screened 


log 2 (n r) 


r l/3 



TABLE I. Large-r behavior of C(r) (C(r) = - log G e (r), 
where G e (r) is the electron Green's function) for the striped 
state and isotropic state of a 2DEG with a half-filled up- 
per Landau level (the results for the isotropic case are from 
Refs. 21 and 22). For both unscreened and screened Coulomb 
interaction, tunneling is weaker in the isotropic case. 





v* w 1/2, L > 1 
striped 


v w 1/2, L = 
isotropic 


unscreened 


exp ( — const / Vv ) 


exp (— const /V) 


screened 


exp(— const log 2 (V)) 


exp ( — const / \/V) 



TABLE II. General behavior of the I-V curves for the 
striped state and isotropic state of a 2DEG with a half-filled 
upper Landau level (the results for the isotropic case are from 
Refs. 21 and 22) . 
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